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PREFACE 

The  study  reported  herein  was  authorized  by  the  Office,  Chief  of 
Engineers  (OCE),  on  22  May  19T5,  at  the  request  of  the  U.  S.  Army  Engi- 
neer Division,  New  England  (NED). 

The  investigation  was  conducted  during  the  period  July  1975  to 
June  1976  in  the  Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  the  direction  of  Messrs.  H.  B. 
Simmons,  Chief  of  the  Hydraulics  Laboratory,  and  J.  L.  Grace,  Jr.,  Chief 
of  the  Structures  Division.  The  study  was  conducted  by  Messrs.  M.  S. 
Dortch,  B.  Loftis,  D.  G.  Fontane,  and  S.  C.  Wilhelms  with  as'sistance 
from  J.  H.  Riley.  This  report  was  prepared  by  Messrs.  Dortch,  Fontane, 
Loftis,  and  Wilhelms  and  reviewed  by  Mr.  Grace. 


Representatives  of  WES  met  with  Mr.  R.  J.  DiBuono  of  NED  at  WES 
diiring  November  1975  and  with  Messrs.  D.  P.  Buelow  and  R.  J.  DiBuono  of 
NED  and  Messrs.  E.  E.  Eiker  and  S.  B.  Powell  of  OCE  during  May  1976  to 
discuss  the  scope,  objectives,  and  approach  of  the  study.  During  the 
course  of  this  study,  Mr.  Charles  Wener  of  NED  was  detailed  to  WES  to 
assist  in  physical  model  testing.  Mr.  D.  P.  Buelow  of  NED  also 
visited  WES  during  June  1976  to  review  the  progress  to  date  and  to  gain 
an  understanding  of  the  mathematical  model.  Messrs.  Fontane  and  Loftis 
visited  NED  during  July  1976  to  deliver  the  mathematical  model  and 
assist  in  preparing  the  model  for  use  in  NED's  computer  facilities. 

Directors  of  WES  dioring  this  study  and  the  preparation  and  publi- 
cation of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L.  Cannon, 
CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 
. UNITS  OF  MEASUREMENT 


U.  S.  customary  xuiits  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Sl)  \mits  as  follows: 


Multiply 

inches 

feet 

miles  (U.  S.  statute) 

square  feet 

acres 

acre-feet 

feet  per  second 

cubic  feet  per  second 

feet  per  second  per  second 

Btu'  s 

Fahrenheit  degrees 




25.1+ 

0.301+8 

1.60931+4 

0.09290304 

4046.856 

1233.482 

0.3048 

0.02831685 

0.3048 

1055.056 
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To  Obtain 

millimetres 

metres 

kilometres 

square  metres 

square  metres 

cubic  metres 

metres  per  second 

cubic  metres  per  second 

metres  per  second  per 
second 

Joules 

Celsius  degrees  or  Kelvins* 


1 


* To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  following  formula:  C = (5/9)(F  - 32).  To  obtain  Kelvin 
(K)  readings,  use:  K = (5/9)(F  - 32)  + 273. 15. 


DICKEY-LINCOLN  SCHOOL  LAKES  HYDROTHERMAL  MODEL  STUDY 


Hydraulic  Laboratory  Investigation 

PART  I : INTRODUCTION 

Purpose  and  Scope  of  Study 

1.  The  overall  objective  of  this  study  was  to  provide  the 

U.  S.  Army  Engineer  Division,  New  England  (NED),  with  an  understanding 
of  the  hydrothermal  processes  unique  to  the  Dickey-Lincoln  School 
Lakes  project  and  to  incorporate  these  processes  into  a numerical  simu- 
lation model  that  NED  could  use  to  predict  temperatures  within  and 
downstream  from  the  lakes  for  planning  and  design  purposes.  Initially, 
it  was  planned  that  the  capability  to  simulate  the  dissolved  oxygen 

(D.O. ) regimes  within  and  downstream  from  the  lakes  by  means  used 
1-h 

in  previous  studies  would  be  included  in  the  numerical  simulation 
technique.  However,  because  of  the  lack  of  adequate  data  relative 
to  estimates  of  the  total  daily  oxygen  depletion  rate  in  large  lakes 
in  this  locality  that  is  needed  for  application  of  the  referenced 
technique  and  the  lack  of  an  adequate  methodology  for  estimating 
1 the  D.O.  uptake  in  the  reregulating  pool,  NED  decided  to  use  exist- 

ing limnological  analysis  techniques  and  the  limited  data  available 
for  estimating  far-fi‘=‘ld  D.O.  profiles  required  relative  to  fishery 
,j  interests. 

2.  Because  of  the  highly  dynamic  nature  of  operational  char- 

, acteristics  of  the  proposed  pumped-storage  project,  physical  hydrau- 

lic models  were  used  to  determine  the  hydrodynamic  response  to 
various  modes  of  operating  the  pimped-storage  hydropower  project. 
Physical  model  hydrodynamic  input  was  necessary  to  develop  a descrip- 
tive mathematical  model.  The  mathematical  model  provided  the  capa- 
bility for  assessing  the  effect  of  historical  data  on  the  lakes  for 
yeeirlong  periods. 
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Pro.lect  Description 
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3.  The  proposed  Dickey-Lincoln  School  hydropower  project  would  be 
located  on  the  St.  John  River  in  northern  Maine.  The  two  damsites, 

Dickey  and  Lincoln  School,  are  approximately  30  miles*  west  of  Fort  Kent, 
Maine.  Dickey  damsite  is  just  upstream  of  the  confluence  of  the  Allagash 
and  St.  John  Rivers  (Figure  l).  Dickey  Lake  would  extend  upstream  about 
miles,  reaching  into  Quebec  Province,  Canada.  Lincoln  School  damsite 
is  about  11  miles  downstream  of  the  Dickey  damsite.  The  lake  would  ex- 
tend to  the  tailwater  of  the  Dickey  Dam  and  for  a,  distance  of  2 miles 
up  the  Allagash  River.  Both  lakes  would  provide  hydropower  d\iring  peak 
power  demand  periods.  Pumped  flow  from  Lincoln  School  Lake  to  Dickey 
Lake  dviring  low  energy  demand  periods  is  proposed  for  the  project  to 
maintain  the  Dickey  pool,  thus  increasing  the  power  production  capabil- 
ities. The  proposed  pumpback  would  pass  through  the  multiple  penstock 
intake  structure. 

k.  Dickey  Dam  would  include  two  earth-fill  embankments,  termed 
the  North  and  South  Dams,  separated  by  a rock  hill  on  the  right  abutment 
of  the  St.  John  River.  The  outlet  works  and  hydroelectric  power  plant 
would  be  located  in  this  hill  (Plate  l).  The  hydraulic  features  of 
Dickey  Dam  consist  of: 

a.  Multipurpose  intake  structure,  with  selector  gates  to 
provide  water-quality  control,  which  discharges  thro\xgh  a 
30-ft-diam  bifurcated  conduit  connected  to  one  turbine  in 
the  powerhouse  and  to  the  flood  control  outlet  works. 

b.  Multiple  penstock  intake  structure,  with  selector  gates, 
which  discharges  through  three  (authorized)  to  five 
(proposed)  27-ft-diam  penstocks  each  connected  to  one  unit 
in  the  powerhouse. 

c_.  Powerhouse  with  fo\ir  conventional  (authorized)  or  three 
conventional  and  one  reversible  (proposed  initial  devtiop- 
ment)  turbines  with  provisions  for  two  additional  re- 
versible turbines  (proposed  ultimate  development). 

d.  A converging  chute  spillway  with  crest  el  910.** 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  presented  on  page  3. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level  datum. 
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Location 


With  the  pool  at  the  maximum  operating  level  (el  910),  the  total  storage 
of  Dickey  Lake  would  be  7,700,000  acre-feet  and  the  surface  area  86,000 
acres.  The  streambed  elevation  at  the  dams it e is  el  580. 

5.  Lincoln  School  Dam  would  consist  of  an  earth-fill  embankment, 
a gated  spillway,  and  powerhouse  with  three  turbines.  The  Lincoln 
School  reregulation  reservoir  would  have  a total  storage  of  about  85,000 
acre-feet  at  a maximum  operating  level  of  el  620. 

Approach 

6.  The  thermal  characteristics  of  lakes  are  affected  by  flow 
processes  such  as  inflow  mixing  and  displacement,  outflow  withdrawal 
distribution,  diffusion,  and  internal  currents.  With  a lake  having  the 
highly  dynamic  flow  conditions  and  flow  magnitudes  associated  with  the 
Dickey-Lincoln  School  project,  an  understanding  of  these  flow  processes 
is  necessary  if  a reasonable  description  of  the  thermal  characteristics 
is  to  be  made. 

7.  Two  physical  hydraulic  models  were  used  to  provide  an  im- 
proved understanding  of  the  hydrodynamics  unique  to  the  project.  An 
undistorted,  1: 200-scale  model  of  the  Dickey  Lake  intake  structures  and 
local  topography  was  used  to  determine  the  steady-state  selective  with- 
drawal and  pumpback  characteristics  for  various  conditions.  A highly 
distorted-scale  model  (l:3600  horizontal,  l:l80  vertical)  was  used  to 
simulate  the  entire  dual  reservoir  system  and  determine  the  response  to 
dynamic,  unsteady-state,  density-stratified  flow.  Information  from  the 
two  physical  models  was  used  to  modify  existing  and  develop  new  al- 
gorithms for  a one-dimensional  mathematical  model. 

8.  The  mathematical  model  allowed  simulation  of  the  hydrodynamic 
and  heat  exchange  characteristics  so  that  the  thermal  regimes  within 
and  downstream  of  the  two  lakes  could  be  evaluated  for  various  hydro- 
logic  and  meteorologic  conditions  and  various  pumped-storage  hydro- 
power  operations. 
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PART  II:  DICKEY  LAKE  INTAKE  MODEL 

Purpose  and  Scope 

9-  The  undistorted,  1: 200-scale  physical  model  was  required  to 
determine  the  steady-state  selective  withdrawal  and  pumpback  mixing 
characteristics  unique  to  the  intake  structures  and  local  topography  of 
Dickey  Lake.  Because  of  the  high  discharge  rates  (i+0, 000-60, 000  cfs) 
and  unusual  local  topography,  a physical  model  was  needed  to  determine 
whether  or  not  the  vertical  flow  distribution  and  outflow  temperature 
could  be  predicted  from  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  Generalized  Selective  Withdrawal  Technique^  and  to  provide 
information  with  which  to  modify  the  predictive  technique  if  discrepan- 
cies were  found.  The  generalized  technique  was  used  to  route  outflow  in 
the  numerical  model.  The  undistorted  model  was  also  needed  to  determine 
the  gross  entrainment,  mixing,  and  dilution  characteristics  of  the  pump- 
back  jet  for  various  modes  of  operation.  This  information  was  required 
for  mathematical  model  improvements. 

Scale  Relations 


10.  The  predominant  forces  affecting  density-stratified  flows  in 
lakes  are  inertia  and  gravity  as  modified  by  density  differences.  In 
such  cases,  hydraulic  similarity  between  a model  and  prototype  system 
requires  that  the  ratio  of  inertial  to  gravitational  forces,  defined  as 
the  Froude  number  of  flow,  be  the  same  in  both  the  model  and  the  proto- 
type. With  the  density  differences  in  the  model  set  equal  to  those  in 
the  prototype,  the  accepted  equations  of  hydraulic  similitude,  based  on 
the  Freudian  relations,  were  used  to  express  the  mathematical  relations 
between  the  dimensional  and  hydraulic  quantities  of  the  model  and  the 
prototype.  The  general  relations  for  transfer  of  model  data  to  proto- 
type equivalents  are  as  follows; 
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Scale 


Dimension 

Ratio 

Relation 

Length 

L 

r 

= 

1:200 

Velocity 

V 

r 

= 

r 

l:li|.lU 

Time 

T 

r 

= 

l1/2 

r 

l:lU.ll+ 

Volume  flow  rate 

% 

= 

l5/2 

r 

1:565,685 

Density  difference 

Ap 

v» 

= 

1 

1:1 

Description 

11.  The  Dickey  Lake  intake  model  (Figure  2)  reproduced  the  pre- 
liminary details  of  the  multiple  penstock  and  multipurpose  intake  struc- 
tures (Plates  2 and  3)  and  UOOO  by  UOOO  ft  of  approach  topography.  The 
structural  members,  piers,  wet  well,  and  bell-mouthed  intakes  of  the 
Dickey  intake  structures  were  constructed  of  transparent  plastic.  The 
selector  gates,  which  were  made  of  sheet  metal,  could  be  raised  and 
lowered  manually.  The  trashracks  were  simulated  with  1/8-in.  wire  mesh. 
The  model  of  the  multiple  penstock  intake  structure  included  five  pen- 
stocks so  the  authorized  (three  27-ft-diam  penstocks)  and  proposed  (five 
27-ft-diam  penstocks)  plans  could  be  tested.  Urethane  foam  was  con- 
toured to  simulate  topography  surrounding  the  structures.  The  model  was 
contained  in  a 20-  by  20-ft  transparent  plastic  flume  so  vertical  ve- 
locity profiles  could  be  easily  observed. 

12.  Saline  and  fresh  waters  were  used  to  reproduce  the  density 
variations  that  are  anticipated^  in  the  prototype  due  to  temperature 
differences.  Density  stratification  was  created  by  placing  fresh  water 
over  saline  water  by  means  of  an  overflow  weir,  the  crest  elevation  of 
which  could  be  varied  to  achieve  the  desired  density  distribution. 
Density  profiles  were  determined  from  temperature  and  conductivity  data 
obtained  with  sensors  that  traversed  the  vertical  direction,  ‘Generation 
and  pumpback  flow  rates  were  regulated  with  hand-operated  valves  and 
measured  with  rotameters.  Pumpback  water  was  pumped  from  a supply  tank 
in  which  a brine  solution  was  prepared.  Vertical  velocity  distributions 
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were  obtained  by  dropping  dye  particles  into  the  fliune  and  recording 
the  displacement  of  the  resulting  dye  streaks  with  video  equipment.  This 
is  a very  reliable  technique  for  determining  velocities  with  magnitudes 
as  low  as  those  occurring  in  lake  currents.  Water-surface  elevations 
were  measured  with  a staff  gage. 


Model  Tests  and  Results 


Selective  withdrawal  tests 

13.  The  anticipated  selective  withdrawal  plan  during  periods  of 

stratification  consisted  of  maintaining  a submerged-weir  flow  condition 

by  withdrawing  water  over  the  top  of  the  selector  gates  (Figure  3). 

This  mode  of  operation  was  tested  with  both  intake  structures  for 

* 

various  stratification  conditions  and  heads  on  the  weirs  H that 

w 

ranged  from  23  ft  to  43  ft.  A 23-ft  head  was  the  minimum  head  allowed 
by  the  project  design.  Control  valves  allowed  the  operation  of  single 
or  multiple  units  of  the  multiple  penstock  intake  structure.  Discharges 
were  set  at  10,^00  cfs  per  \init;  releases  from  all  six  units  totaled 
60,000  cfs. 

14.  Density  profiles  were  obtained  before  each  test.  Velocity 
profiles  and  outflow  density  were  recorded  during  the  test.  By  providing 
the  density  stratification,  weir  crest  elevation,  weir  length,  discharge 
rate,  and  pool  elevation,  selective  withdrawal  predictions  were  made 
with  the  WES  genera] ized  technique.  The  generalized  predictive  technique 
has  been  coded  for  computer  application  and  entitled  SELECT.  A version 
of  SELECT  is  also  included  as  an  outflow  subroutine  to  the  mmerical 
simulation  model  described  in  Part  IV.  Each  prediction  was  compared 
with  the  observed  velocity  distribution  and  outflow  density  to  verify 
the  accuracy  or  determine  the  inaccuracy  of  SELECT.  Where  discrepancies 
existed  between  physical  model  and  generalized  results,  physical  model 
results  were  used  to  modify  portions  of  SELECT  to  provide  a more 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  B). 
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Definition  sketch  of  selective  withdrawal 
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I acclirate  description  of  the  withdrawal  characteristics  unique  to  this 

t ■ 

i project. 

j Selective  withdrawal  results 

^ 15.  Releases  from  the  multipurpose  intake  structure  and  any  one 

i 

^ imit  of  the  miiltiple  penstock  intake  structure  produced  similar  with- 

[.  drawal  characteristics.  These  characteristics  correlated  well  with 

SELECT  predictions  for  conditions  of  mild  stratification  (spring  or  fall) 
or  greater  than  40  ft.  However,  for  summer  stratification  condi- 

tions and  H less  than  40  ft,  a greater  contribution  of  flow  from  the 
w 

lower  portion  of  the  velocity  profile  and  a larger  outflow  density  (about 
3°F  colder  for  a 30°F  temperature  difference  from  surface  to  bottom)  were 
observed  than  that  predicted  by  SELECT.  The  increased  contribution  of 
flow  from  below  the  location  of  the  maximum  velocity,  V , was 

C18.X 

attributed  to  a slightly  superelevated  thermocline  immediately  upstream 
of  the  weirs  resulting  from  a hydrostatic  pressure  gradient  reduction 
near  the  weir  crest  caused  by  local  acceleration.  With  the  weir  for 
each  wet  well  segmented  into  four  individual  gates  separated  by  piers, 
contraction  of  the  flow  around  the  piers  could  cause  local  accelerations. 
Also,  high  rates  of  discharge  (10,000  cfs  per  penstock)  over  the  weirs, 
relatively  small  with  respect  to  the  discharge,  and  the  position 

of  the  thermocline  near  the  weir  crest  could  contribute  to  these  accel- 
erations. These  accelerations  would  certainly  be  more  prono\uiced  for 
small  and  strong  stratification.  This  is  a reasonable  explanation 

for  the  fairly  small  discrepancy  between  the  model  and  generalized  re- 
sults^ (select)  considering  that  the  generalized  results  were  obtained 
for  a uniform,  continuous,  single  weir  that  extended  the  full  width  of 
a 1-ft-wide  flume. 

l6.  A distribution  of  flow  similar  to  that  discussed  above  was 
observed  for  simultaneous  releases  through  three  or  more  units  of  the 
multiple  penstock  intake  structure  for  the  entire  range  of  stratifica- 
tion and  H conditions  tested.  Additionally,  the  elevation  of  V 

w xnsoc 

and  th^  elevation  of  the  lower  limit  of  withdrawal,  , were  observed 

to  be  lower  than  the  predicted  elevations  for  releases  through  three  or 
more  of  these  five  penstocks.  These  deviations  in  the  ant-icipated  flow 
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J phenomena  were  believed  to  be  created  by  topographical  approach  con- 

ditions. It  was  concluded  that  the  approach  geometry  and  elevated 
topography  (knoll)  combined  with  high  discharges  created  a restricted 
(converging)  flow  field  upstream  of  the  weirs.  This  flow  restriction  in- 
fluenced the  withdrawal  characteristics  farther  out  in  the  lake.  Al- 
^ though  the  approach  topography  of  the  multiple  penstock  intake  structure 

did  not  completely  control  the  withdrawal  characteristics,  it  did  have 
an  increasing  influence  as  the  number  of  units  operating  was  increased 
from  three  to  five.  With  all  five  units  of  the  multiple  penstock  intake 
structure  operating,  the  maximum  variation  in  observed  temperatures  (pro- 
jected from  densities)  of  outflow  from  predicted  values  was  about  6°F 
colder  or  20  percent  of  the  total  temperature  difference  expected  in  the 
pool  during  the  summer.  The  influence  of  the  topography  was  negligible 
when  operating  a single  unit  of  the  multiple  penstock  structure,  and  as 
stated  previously,  the  withdrawal  characteristics  compared  closely  with 
those  of  the  multipurpose  intake. 

Selective  withdrawal  modifi- 
cations  for  mathematical  model 

IT • Three  basic  revisions  were  made  to  the  mathematical  model  to 
improve  the  accuracy  of  withdrawal  predictions.  These  revisions  con- 
sisted of:  (a)  reduction  of  the  actual  weir  length  to  force  a lower 

limit  of  withdrawal,  (b)  lowering  of  the  maximum  velocity  elevations, 
and  (c)  modification  of  the  equation  that  predicts  the  velocity  distri- 
bution below  V . The  specifics  of  each  revision  are  discussed  in 
max 

the  following  paragraphs. 

l8.  To  account  for  the  deepened  withdrawal  zone  created  by  the 
topographical  flow  restriction  upstream  of  the  multiple  penstock  struc- 
ture during  multiunit  operation,  modifications  were  made  to  the  mathe- 
matical model  to  reduce  the  weir  length,  thus  decreasing  the  elevation 

of  the  lower  limit.  The  lower  withdrawal  limit,  Z , is  calculated 

o 

t from  the  equation 


ll* 


(1) 


where 


) 


V = average  velocity  over  the  weir,  fps 
w 

Z = vertical  distance  from  the  elevation  of  the  weir  crest  to 
the  lower  limit  of  the  zone  of  withdrawal,  ft 

H = head  on  the  weir  for  free  flow  or  depth  of  flow  over  the 
weir  for  submerged  flow,  ft 

Ap  = density  difference  of  fluid  between  the  elevations  of  the 

weir  crest  and  the  lower  limit  of  the  zone  of  withdrawal, 
g/cc 

p = density  of  fluid  at  the  elevation  of  the  weir  crest,  g/cc 

2 

g = acceleration  due  to  gravity,  ft/sec 


Because  V is  calculated  from 
w 


V 

w 


bH 

w 


(2) 


where 

Q = discharge  over  the  weir,  cfs 
b = weir  length 

a reduction  in  weir  length  causes  an  increase  in  Z^  . The  effective 
weir  length  corresponding  to  the  number  of  units  operating  is  provided 
in  the  following  tabulation: 


Number  of 

Actual  Weir 

Effective  We 

Units  Operating 

Length,  ft 

Length,  ft 

1 

80 

80 

2 

160 

160 

3 

2i*0 

230 

1+ 

320 

290 

5 

i+00 

3i+0 

6 

U80 

380 

For  the  majority  of  the  tests  conducted,  the  elevation  of  the  lower 
limit  ranged  from  el  700  to  800. 

19.  The  equation  used  in  SELECT  to  calculate  the  height  of  the 
maximum  velocity  is 


H + Z 


sin  1.57 


H + Z 
w c 


t2 


(3) 
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where 

= vertical  distance  from  the  elevation  of  the  maximum  velocity 
to  the  lower  limit  of  withdrawal,  ft 

With  three  or  more  units  of  the  multiple  intake  structure  in  operation, 
the  location  of  the  maximum  velocity  was  observed  to  be  lower  in  the 
pool  than  was  predicted  by  Equation  3.  This  phenomenon  was  not  observed 
for  higher  heads  (H^  greater  than  33  ft).  With  of  23  ft,  the  lo- 

cation of  Y^  was  found  to  be  about  10  ft  lower  than  predicted.  The 
mathematical  model  was  modified  to  calculate  a revised  location  of  the 
maximum  velocity  with  respect  to  the  lower  limit  of  withdrawal  from 


where  Y^  is  calculated  from  Equation  3.  Equation  U is  only  used  if 
three  or  more  units  of  the  multiunit  structure  are  operating  and 
is  less  than  33  ft. 

20.  The  equation  obtained  from  the  WES  generalized  selective 
withdrawal  results^  for  describing  the  dimensionless  velocity  distri- 
bution for  the  portion  below  the  maximum  velocity  with  submerged  weir 
flow  is 


where 


\"Plm 


v^  = local  velocity  in  the  zone  of  withdrawal  at  a distance  y 

below  the  elevation  of  the  maximum  velocity  V , fps'^ 

max  ’ ^ 

y^  = vertical  distance  from  the  elevation  of  the  maximum  veloc- 
ity Vmax  to  that  of  the  corresponding  local  velocity 
v^  , ft 

V = maximum  velocity  in  the  zone  of  withdrawal,  fps 
max  ’ ^ 

Ap^  = density  difference  of  fluid  between  the  elevations  of  the 

maximum  velocity  and  the  corresponding  local  velocity 


. g/cc 


Ap^  = density  difference  of  fluid  between  the  elevations  of  the 

maximum  velocity  and  the  lower  limit  of  the  zone  of  with- 
drawal, g/cc 
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Equation  5 was  modified  as  follows  to  accoimt  for  situations  in  which  a 
larger  contribution  of  flow  was  realized: 


V^lm/ 


(6) 


where  p varied  from  1.8  to  3.0. 

Equation  6 was  incorporated  into  the  mathematical  model  to  be  used  as 
follows : 

Whenever  three  or  more  units  of  the  multiple  intake 
structure  are  operating,  p = 1.8  for  all  conditions. 

b.  For  fewer  than  three  of  the  multiple  intakes  discharging, 
during  summer  stratification,  ana  H less  than  40  ft, 
p is  calculated  from 

p = 1.8  + 0.0T06(H  - 23.0)  (T) 

w 

For  the  same  conditions,  except  with  H of  UO  ft  or 
greater,  p = 3.0  . 

c_.  For  fewer  than  three  of  the  multiple  intakes  dis- 
charging d’oring  mild  stratification  (spring  or  fall), 

P = 3.0  . 

21.  The  exponent,  p , of  1.8  was  obtained  by  fitting  model  data 
(obtained  from  test  conditions  of  single  unit  operation  with  siimmer 
stratification  and  H = 23  ft  and  condition  a above)  to  the  form  of 

W “ 

Equation  6,  using  the  method  of  least  squares.  The  data  fit  gave  a 
correlation  coefficient  of  0.95  and  a standard  error  of  0.08  for  predic- 
tions of  V /V  . These  model  data,  the  least-squares  curve  fit 
X ms>x 

(p  = 1.8),  and  Equation  6 with  p = 3.0  (WES  generalized  results)  are 
all  plotted  in  Plate  h for  comparison  purposes. 

Pumpback  tests 

22.  Basically,  two  modes  of  operation  were  tested  for  pumpback 

characteristics:  pumping  over  the  top  of  the  selector  gates  (weirs) 

and  pumping  under  the  gates  with  the  gates  in  the  raised  position 
(Figure  U).  Variations  in  the  pumpback  density,  the  density 
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stratification,  and  the  number  of  units  operating  were  investigated  for 
both  configiirations. 

23.  The  layers  of  the  lake  that  contributed  to  the  entrainment 
process  were  determined  by  visual  observation,  and  the  amount  of  entrain- 
ment was  determined  from  density  measurements  and  dilution  calculations. 
A definition  sketch  of  the  pumpback  characteristics  is  provided  in  Fig- 
ure 5*  For  each  test  condition,  an  initial  density  profile  (prior  to 
pumpback),  the  density  of  the  pumped  water,  , the  average  density 

of  the  pumpback  current,  , and  the  final  density  profile  (after  ter- 

mination of  pumped  flow)  were  measured.  The  initial  and  final  density 
profiles  were  plotted  to  determine  the  approximate  average  density  of 
the  pumpback  current,  p^  . An  example  of  this  technique  is  shown  in 
Figure  6.  The  p^  value  was  compared  with  measurements  taken  in  the 
density  current  d\iring  pumpback  to  substantiate  the  procedure.  These 
two  densities  usually  correlated  very  well. 

2k,  Conservation  of  mass  and  volume  was  used  to  determine  the 
amount  of  gross  entrainment  as  follows.  From  the  conservation  of  mass 


pQ  = pQ,  + pQ 
c c 00  e e 


where 

p^  = average  density  of  the  pumpback  current,  g/cc 
= volume  flow  rate  of  the  pumpback  cirrrent,  cfs 
p^  = density  of  the  pumped  water,  g/cc 

= volume  flow  rate  of  water  pumped,  cfs 
p^  = average  density  of  the  entrained  flow,  g/cc 
= volume  rate  of  entrained  flow,  cfs 
Also,  from  the  continuity  of  volume  for  an  incompressible  fluid 


% - ^o  ^e 


Substitution  of  Equation  9 into  Equation  8 for  yields 


^e  = ^«o 
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Figure  5.  Definition  sketch  of  pumpback  characteristics 


Figure  6,  Effect  of  pumpback  on  density 
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where 


e = the  entrainment  coefficient 


and 


e = 


_o c_ 

P - P 
c e 


The  value  of  was  taken  from  the  initial  density  profile  from  the 

elevation  at  which  the  maximum  contribution  of  entrainment  was  observed 
in  the  model. 

Pumpback  results 

25.  The  gross  entrainment  coefficient,  e , was  determined  for 
each  pumpback  test  condition.  The  average  values  of  e for  the  various 
configurations  tested  are  tabulated  below. 

Pumping  Buoyancy 


Configuration 

of  Jet 

e 

Over  the  gates 

Negative 

0.3 

Over  the  gates 

Positive 

0.3 

Under  the  gates 

Positive 

0.6 

With  pumpback  over  the  gates  and  negative  buoyancy,  water  was  entrained 
into  the  jet  primarily  from  the  layers  of  the  pool  in  the  vicinity  of 
and  above  the  elevation  of  the  weir  crest.  With  pumpback  over  the  gates 
and  positive  buoyancy,  water  was  entrained  into  the  jet  primarily  from 
the  layers  of  the  pool  existing  below  the  weir  crest  down  to  el  8OO. 

The  maximum  contribution  was  observed  to  be  about  20  ft  below  the  weir 
crest.  With  pumpback  under  the  gates  (buoyancy  in  this  case  was  always 
positive  because  of  the  depth  of  the  penstock  center  line,  el  823.5), 
water  was  primarily  entrained  from  layers  between  el  8OO  and  85O.  The 
higher  e for  this  configuration  was  attributed  to  the  greater  veloc- 
ity of  the  jet  exiting  from  the  bell  mouth  of  the  penstock.  For  pump- 
back over  the  weir,  the  flow  was  able  to  spread  out  to  the  full  weir 
length,  causing  a smaller  average  velocity  of  the  jet  and  a smaller  e . 
The  model  indicated  that  the  number  of  units  pumping  back  had  a negli- 
gible effect  on  e , probably  because  the  unit  discharge  and  average  ve- 
locity of  the  jet  were  the  same  for  three  units  pvunping  as  for  one. 
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Pumpback  algorithm 
for  mathematical  model 

26.  The  mathematical  model  was  modified  to  include  the  pumpback 
results.  To  facilitate  the  description  of  entrainment  in  the  mathemati- 
cal model,  entrainment  was  assumed  to  occur  from  only  one  layer  of  the 
pool.  The  location  of  this  layer  is  dependent  on  the  pumping  configura- 
tion and  buoyancy  of  the  pumpback  jet.  For  pumpback  over  the  gates  with 
a negative  buoyancy  (p^  greater  than  density  existing  in  the  pool  at 
the  elevation  of  the  weir  crest,  P^  'the  layer  that  includes 

p are  used  for  entrainment.  For  pumpback  over  the  gates  with  a posi- 
w 

tive  buoyancy  (p  < p ),  water  is  entrained  from  the  layer  20  ft  below 
o w 

the  elevation  of  the  weir  crest.  For  pumpback  under  the  gates,  water  is 
entrained  from  the  layer  that  includes  el  823,  which  is  about  the  center 
line  of  the  penstock.  Water  from  the  entrainment  layer  is  mixed  with 
the  pumped  discharge  in  the  amount  of  eQ^At  , where  At  is  the  pumping 
period.  The  density  of  this  mixture,  p^  , is  determined  and  the  pump- 
back mixture  is  placed  into  the  lake  at  the  elevation  that  p^  exists 
in  the  density  profile.  Further  descriptions  of  budgeting  and  mixing 
related  to  pumpback  are  discussed  in  the  hydrodynamic  and  mathematical 
model  sections  (Parts  III  and  IV,  respectively). 
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PART  III:  LAKE  HYDRODYNAMICS  MODEL 


Purpose 

27.  The  purpose  of  the  distorted-scale  model  was  to  aid  in  de- 
fining the  hydrodynamics  of  Dickey-Lincoln  School  Lakes  resulting  from 
the  expected  operating  conditions  of  the  proposed  prototype.  The  model 
was  needed  to  determine  the  effects  of  unsteady  generation  and  pumpback 
on  the  density  stratification  of  Dickey  and  Lincoln  School  Lakes  during 
any  particular  stratification  period  and  to  help  identify,  validate, 
and  quantify  any  modifications  to  be  made  to  the  mathematical  model  to 
improve  the  reliability  of  the  predictions. 

Scale  Relations 


28.  The  accepted  equations  of  hydraulic  similitude,  also  based  on 
the  Froudian  criteria,  were  used  to  express  the  mathematical  relations 
between  dimensions  and  hydraulic  quantities  of  the  model  and  prototype. 
Allowing  for  vertical  scale  distortion,  the  general  relations  for  trans- 
fer of  model  data  to  prototype  equivalents  are  as  follows: 


Dimension 

Ratio 

Scale  Relation 

Length  in  vertical  direction 

L 

r 

= 

L 

y 

I:l80 

Length  in  horizontal  direction 

L 

r 

= 

L 

X 

1 : 3600 

Area  in  a vertical  plane 

A 

L L 

0 

0 

0 

CO 

r 

X y 

Area  in  a horizontal  plane 

A 

r 

= 

l2 

X 

1:12,960,000 

Velocity 

V 

r 

= 

L^/2 

y 

1:13. 

Time 

T 

_ 

L /L^/2 

1:268 

r 

X y 

P'low 

Q 

— 

L l3/2 

1:8,693,832 

r 

X y 

Density  difference 

= 

1 

1:1 

Measurements  of  flow,  water-surface  elevations,  and  time  can  be  trans- 
ferred quantitatively  from  the  model  to  the  prototype  by  means  of  the 
scale  relations  above. 


Description 


29.  The  model  (Figure  7)  was  constructed  to  highly  distorted 
length  scale  ratios  of  l;l80  vertically  and  1:3600  horizontally.  With 
the  vertically  distorted  model,  turbulent  flow  was  preserved  while  sim- 
ulating the  entire  dual  reservoir  system  in  order  to  maintain  the  same 
fundamental  character  of  flow  as  in  the  prototype.  Use  of  an  undis- 
torted model  of  the  entire  system  would  require  a model  of  such  large 
dimensions  and  discharge  capacity  that  the  cost  would  be  impractical. 

By  horizontal  scale  compression,  the  total  size  of  the  model  was  reduced 
and  hydraulic  similitude  was  preserved. 

30.  The  model  lakes  were  constructed  of  transparent  and  opaque 
plastic  to  facilitate  photography  and  visual  observations  of  currents 
in  the  vertical  and  longitudinal  directions.  The  model  approximated 
the  geometry  and  reproduced  the  scaled  elevation-storage  relationships 
of  the  prototype  lakes.  The  sides  of  the  model  flvune  were  stepped 
vertically  (Figure  8)  to  satisfy  the  above  requirements  without  visual 
distortion  through  the  side.  Urethane  foam  was  contoured  to  simiilate 
the  knoll  where  the  intake  structures  are  to  be  located,  the  topography 
at  the  confluence  of  the  Allagash  and  St.  John  Rivers,  the  configuration 
of  the  tailrace  of  Dickey  Dam,  and  a large  island  located  in  the  Little 
Black  arm  of  Dickey  Lake. 

31.  As  in  the  undistorted  model,  the  anticipated  density  strat- 
ification was  generated  with  saline  and  fresh  waters,  and  density 
measurements  were  performed  with  conductivity  and  temperature  sensors. 
Velocity  distributions  were  again  determined  from  video  recordings  of 
dye  streak  displacement.  Various  dyes  were  used  to  follow  the  movement 
of  water  of  particular  interest. 

32.  Because  of  the  compressed  horizontal  scale,  it  was  not  prac- 
tical to  reproduce  the  2k  gates  and  6 intakes  of  the  intake  structures 
in  the  scaled  horizontal  dimension  of  0.13  ft.  Therefore,  the  selector 
gates  were  combined  into  a single  gate  that  was  vertically  controlled 
by  an  electromechanical  actuator  and  the  six  intakes  were  combined  into 
two.  Provisions  for  two  intakes  and  wet  wells  were  necessary  to  maintain 


Dickey-Lincoln  School  Leikes  hydrodynamic  model 


MODEL  CROSS  section 


the  proper  (scaled)  inlet -outlet  dimensions  used  during  the  two  differ- 
ent flow  processes  of  generation  and  pumphack.  Generation  and  pumpback 
hydrographs  of  Dickey  Lake  were  simulated  by  means  of  a programmable, 
reversible  flow,  rotary  pump.  Allagash  River  inflow  was  provided  and 
regulated  by  a programmable  rotary  pump.  Inflow  to  Dickey  Lake  was  sup' 
plied  by  a head  tank  and  regulated  with  a rotameter  and  hand  valve. 
Generation  periods  at  Lincoln  School  Dam  were  simulated  by  means  of  an 
automatically  activated  solenoid  valve  which  allowed  water  to  be  re- 
leased from  Lincoln  School  Lake  through  a gate  valve,  which  was  preset 
with  a rotameter  to  give  the  desired  discharge. 

33.  Model  control  devices  could  be  operated  manually  or  auto- 
matically from  the  control  panels  (Figure  9).  Automatic  control  was 
achieved  through  analog  signals  generated  by  magnetic  card  readers  that 
tracked  input  data.  Automatic  control  allowed  unsteady  operation  to 
proceed  for  extended  simulation  periods  (L  to  6 weeks  prototype  time). 
Since  the  physical  model  did  not  have  the  capability  to  reproduce  sur- 
face heat  exchange,  realistic  prototype  simulations  exceeding  6 weeks 
were  not  attempted  because  the  meteorological  effects  could  exceed  the 
hydrodynamic  effects. 
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Model  Tests  and  Results 

Description  of  tests 

31*.  Anticipated  spring,  summer,  and  fall  thei-mal  stratification 
conditions^  (Figure  10)  of  Dickey  Lake,  without  power  operations,  were 
tested  in  the  model.  The  spring  and  fall  conditions  represented  a 5^°F 
water  temperature  in  the  epilimnion  and  h0°¥  in  the  hypolimnion,  and  the 
summer  condition  represented  a 70°F  epilimnion  and  a Uo°F  hypolimnion. 
With  . le  Dickey  pool  at  el  907,  the  thermocline  was  set  at  approximately 
el  900  and  89O  for  the  spring  and  summer  conditions,  respectively.  For 
the  fall  condition,  the  pool  was  set  at  el  903  and  the  thermocline  at 
about  el  860.  Because  the  physical  hydrodynamic  model  did  not  reproduce 
surface  heat  exchange,  the  period  of  time  simulated  by  the  model  was 
relatively  short  so  that  the  hydrodynamic  effects  would  have  a much 
greater  influence  on  stratification  changes  than  would  the  surface  heat 
exchange.  Time  simulated  by  the  model  was  limited  to  2 prototype  weeks 
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Figure  10.  Stratification  conditions  simulated 
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for  spring  and  fall  stratification  and  U to  6 weeks  for  summer  stratifi- 
cation. The  longer  simulation  period  of  the  summer  conditions  was  per- 
missible because  the  rate  of  change  in  density  stratification  due  to 
surface  heat  transfer  is  less  in  the  summer  than  in  the  spring  or  fall. 

35.  The  Lincoln  School  Lake  model  was  stratified  with  fresh  and 
saline  water  to  simulate  a summer  condition  of  75°F  water  for  the  epi- 
limnion  and  65°F  for  the  hypolimnion.  This  was  determined  through  math- 
ematical simulations^  to  be  the  strongest  thermal  stratification  that 
could  be  expected  to  exist  in  Lincoln  School  Lake  with  the  normal 

St.  John  and  Allagash  River  flows.  Tests  were  conducted  to  determine  if 
the  Lincoln  School  water  would  become  fully  mixed  for  the  expected  power 
releases . 

36.  The  model  was  tested,  under  all  stratifications,  on  a steady- 
state  flow  basis  to  evaluate  the  selective  withdrawal  and  pumpback 
characteristics  of  the  intakes.  Similarity  of  these  flow  characteris- 
tics to  the  undistorted  model  results  was  necessary  to  assure  valid 
simulations.  Release  densities,  withdrawal  velocity  profiles,  and  en- 
trainment coefficients  were  determined  for  comparison  with  results  from 
the  undistorted  model  and  SELECT. 

37.  The  effects  of  xuisteady  operation  on  Dickey  Lake  stratifica- 
tion were  investigated  by  operating  the  hydrodynamic  model  on  the  antic- 
ipated prototype  schedules.  Generation  and  pumpback  densities  were 
monitored  to  evaluate  transient  effects  on  release  density  resulting 
from  the  unsteady  operation.  Changes  in  the  stratification  of  Dickey 
Lake  were  evaluated  by  comparing  vertical  density  profiles  taken  at  the 
end  of  each  prototype  week  of  simulation.  The  two  proposed  operation 
plans  tested  are  summarized  in  Table  1.  Both  plans  include  pumpback. 

The  authorized  plan  (no  pumpback)  was  not  tested  because  the  complexity 
of  the  water  mixing  and  transport  characteristics  of  the  project  were 
greatly  reduced,  and  description  of  the  advection  processes  in  the  math- 
ematical model  did  not  require  these  additional  tests. 

38.  During  summer  stratification,  steady-state  inflow  tests  were 
conducted  to  estimate  the  time  of  travel  for  St.  John  River  flows  enter- 
ing Dickey  Lake.  Flow  was  introduced  into  Dickey  headwater  at  an 
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average  rate  of  3200  cfs.  The  time  required  for  the  flow  to  reach  Dickey- 
Dam  and  be  released  was  observed  to  be  90  days  with  the  3200-cfs  steady- 
state  flow. 

39-  A unique  feature  of  the  reregulation  pool,  besides  the  high 
generation  and  pumpback  flow  rates,  was  the  contribution  of  flow  from 
the  Allagash  River.  Because  the  confluence  of  the  Allagash  and  the  re- 
regulation pool  is  in  the  vicinity  of  the  Dickey  tailrace,  the  Allagash 
River  was  considered  to  have  a high  potential  for  contributing  to  the 
quality  of  the  pumpback  flow  from  Lincoln  School  Lake.  Therefore,  tests 
were  conducted  to  determine  what  proportion  of  the  pumpback  water  was 
Allagash  water.  Fluorescent  dye,  sulpha  Rhodamine  B,  -was  used  to  trace 
the  Allagash  water.  With  the  model  simulating  the  unsteady  generation- 
pumpback  modes,  samples  of  pumpback,  generation,  and  Allagash  water  were 
collected  and  analyzed  for  fluorescent  dye  concentration.  The  contribu- 
tion of  the  Allagash  to  the  pumpback  was  determined  from  the  dye  concen- 
trations (paragraph  ) . Allagash  flows  of  1570  and  ^710  cfs  and  un- 
steady operations  described  in  Table  1 were  used  for  these  tests. 
Steady-state  test  results 

UO.  The  steady-state  selective  withdrawal  characteristics  of  the 
distorted-scale  lake  model  compared  closely  with  the  generalized  selec- 
tive withdrawal  results^  but  differed  from  the  undistorted-scale  model 
results.  With  the  6 intakes  and  2k  weirs  combined  into  2 intakes  and  a 
single  weir  in  the  distorted-scale  model,  the  selective  withdrawal 
characteristics  of  the  distorted-  and  undistorted-scale  models  could  be 
expected  to  differ  some.  Also,  because  of  the  distortion,  the  effect 
that  the  approach  topography  had  on  selective  withdrawal  in  the  undis- 
torted model  did  not  exist  in  the  distorted  model.  For  the  simplified 
representation  of  the  intake  structures  and  topography  modeled  in  the 
distorted  model,  the  distorted  model  results  were  expected  to  agree  with 
SELECT  predictions.  Because  the  selective  withdrawal  characteristics 
of  the  distorted-  and  undistorted-scale  models  were  in  reasonable  agree- 
ment, results  from  the  hydrodynamic  model  were  considered  to  be  reliable 
for  understanding  the  dynamic  response  of  the  lakes. 

Ul.  Entrainment  coefficients  obtained  from  the  distorted-scale 
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model  during  steady-state  pumpback  tests  compared  closely  with  those 
obtained  from  the  undistorted-scale  model,  and  adjustment  of  distorted 
model  results  was  not  necessary. 

k2.  With  the  St.  John  River  inflow  into  Dickey  Lake  equal  to  a 
constant  outflow  of  3200  cfs  at  Dickey  Dam  (paragraph  38),  the  dyed  in- 
flow was  observed  to  travel  to  the  dam  in  about  90  prototype  days.  This 
long  travel  time  demonstrated  that  the  effect  of  inflow  on  temperature 
stratification  in  the  vicinity  of  the  dam  was  insignificant  relative  to 
meteorological  effects  in  this  large  lake. 

Results  of  unsteady  operation  and 
mathematical  model  implementation 

)*3.  Releases  from  both  of  the  proposed  hydropower  operation  modes 
{k  and  6 units)  almost  totally  destroyed  the  stratification  setup  in 
Lincoln  School  Lake  (paragraph  35)  within  one  week  of  prototype  opera- 
tion. Only  the  top  2 ft  of  pool  at  the  dam  were  not  well  mixed  after 
one  week  of  prototype  operation.  After  two  weeks  of  prototype  simula- 
tion, the  reregulation  pool  was  totally  mixed.  It  is  reasonable  to  be- 
lieve that  stratification  will  not  develop  in  Lincoln  School  Lake  with 
power  operation.  Therefore,  Lincoln  School  Lake  was  considered  fully 
mixed  for  mathematical  model  application.  Power  releases  from  Dickey 
Lake  were  observed  to  travel  from  the  Dickey  tailrace  to  the  Lincoln 
School  Dam  in  a period  of  approximately  a day.  This  meant  that  net 
flows  into  Lincoln  School  Lake  could  be  mixed  with  water  remaining  in 
the  lake  from  the  previous  day  for  the  mathematical  simulations. 

Water  from  the  Allagash  River  made  a considerable  contribu- 
tion to  pumpback  flow  from  Lincoln  School  Lake  to  Dickey  Lake.  The 
amount  of  Allagash  water  pumped  back  was  found  to  be  a function  of  the 
pumpback  flow  rate  and  the  Allagash  River  Inflow  rate.  With  the  pro- 
posed plan  of  three  reversible  turbines  (pumpback  flow  of  20,000  cfs 
for  6 hr),  the  pumped  discharge  consisted  of  about  20  and  55  percent 
Allagash  River  water  for  Allagash  River  inflows  of  1570  and  U710  cfs, 
respectively.  With  the  alternate  plan  of  one  reversible  tiLrbine  pump- 
ing 7000  cfs  for  6 hr,  the  Allagash  River  water  composed  about  5^  and 
90  percent  of  the  water  pumped  from  Lincoln  School  Lake  for  the  same 
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respective  Allagash  inflow  rates.  This  operational  mode  provided  a 
higher  percentage  of  Allagash  water  than  did  the  3-unit  pumpback  mode, 
probably  because  the  generation  flow  of  i*0,000  cfs  prior  to  pumpback 
with  one  unit  did  not  create  as  high  a backwater  stage  (backing  the 
Allagash  flow  upstream)  as  did  the  generation  flow  of  60,000  cfs.  Also, 
the  ratio  of  the  rate  of  Allagash  River  flow  with  respect  to  the  total 
pumpback  flow  rate  was  higher  for  the  single  pumpback  mode,  thus  in- 
creasing the  influence  of  the  Allagash  River  on  pumpback.  An  equation 
was  developed  to  compute  the  amount  of  Allagash  water  piamped  back  as  a 
percentage  of  the  total  pumped  flow  occurring  in  a pumpback  period. 

This  percentage  of  Allagash  water  to  the  total  pumped  flow,  a , is  a 
function  of  the  Allagash  River  inflow  and  pumpback  flow  rate  and  is 
calculated  from 

0 = (U2.0  - 0.00179Qq)Q^^°'^^  ^o^  (ll) 

where 

= piunpback  flow  rate,  cfs 

= Allagash  River  inflow,  1000  acre- ft /day 
Because  only  two  pumpback  rates  were  tested,  the  application  of  Equa- 
tion 11  should  be  limited  to  those  conditions.  The  Allagash  water  with 
corresponding  temperature  is  mixed  with  either  Lincoln  School  Lake  water 
or  water  released  from  Dickey  Lake  during  the  previous  generation  period 
(if  generation  occurred  on  that  day)  to  obtain  the  temperature  of  water 
pumped  back.  To  conserve  mass  and  energy,  there  is  a limit  to  the  maxi- 
mum amount  of  Allagash  water  that  can  be  pumped  back;  that  is,  pumpback 
water  should  not  exceed  the  total  volume  of  Allagash  flow  that  enters 
Lincoln  School  Lake  in  a day. 

^45.  With  the  distorted-scale  model,  the  pumpback  characteristics 
in  Dickey  Lake  were  observed  for  an  unsteady  operation  plan  over  an  ex- 
tended time  period  (2  to  6 prototype  weeks).  The  flow  conditions  in 
Dickey  Lake  resulting  from  pumpback  with  a submerged  turbulent  jet  con- 
sisted basically  of  two  parts:  a tarbulent  mixing  zone  in  the  vicinity 

of  the  dam  and  gravity-driven  density  currents  in  the  far  field  (Fig- 
ure 11).  The  turbxilent  mixing  zone  extended  from  the  dam  to  a distance 
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Figure  11.  Description  of  pumpback  current 

upstream  of  the  dam  ranging  from  0.5  to  1.5  miles.  Upon  completion  of 
each  pumpback  period,  the  turbulent  zone  would  practically  dissipate 
before  initiation  of  generation  and  form  a mixed  layer  similar  to  the 
convective  current  farther  out  in  the  lake. 

U6.  The  model  demonstrated  the  pumpback  current  to  be  a two- 
dimensional  flow  phenomenon,  which  would  be  expected  for  a long,  narrow 
reservoir.  The  current  was  thickest  at  the  dam  and  gradually  diminished 
in  thickness  from  the  dam  to  the  front  of  the  current,  which  progressed 
up  the  lake.  The  existence  of  a strong  two-dimensional  effect  was  dem- 
onstrated by  the  long  time  required  for  the  current  to  progress  up  the 
lake.  For  summer  conditions  and  with  three  reversible  units  pumping 
back  over  the  weirs,  the  current  traveled  about  half  the  total  distance 
up  the  Little  Black  arm  and  a third  of  the  total  length  of  the  Saint 
John  arm  during  a one-month  simulation.  This  produced  a longitudinal 
variation  in  the  temperature  (density)  profile  within  the  zone  of  inter- 
flow of  the  convective  current.  The  most  significant  aspect  of  this 
longitudinal  variation  is  that  the  density  profile  near  the  dam,  which 
affects  the  outflow  quality  of  Dickey  Lake,  endured  more  mixing  than 
that  farther  out  in  the  lake. 

^^7.  After  completion  of  the  pumpback  period  and  dissipation  of 
turbulence,  the  thickness  and  density  variations  within  the  pumpback 
current  were  fairly  uniform  from  the  dam  to  a distance  of  about  3 miles 
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upstream  of  the  dam.  From  the  physical  model  results,  special  mixing 
procedures  were  developed  for  the  mathematical  model  to  describe  the 
temperature  changes  created  by  the  pumpback  current  in  this  reach  of 
reservoir  adjacent  to  the  dam.  To  maintain  the  energy  and  mass  budget, 
the  advection  process  throughout  the  lake  was  handled  in  a conventional 
manner;  the  pumpback  current  is  placed  in  the  layer  whose  density  most 
nearly  equals  the  density  of  the  pumpback  current. 

The  special  mixing  procedures  discussed  above  consisted  of 
three  parts:  prediction  of  the  pumpback  current  thickness  in  the  vicin- 

ity of  the  dam,  determination  of  mixing  coefficients,  and  partial  mixing 
of  layers  within  the  zone  of  the  current.  For  any  particular  stratifi- 
cation condition  and  operational  mode,  there  was  a corresponding  thick- 
ness of  the  pumpback  current  (zone  of  interflow)  near  the  dam.  The 
thickness  of  this  zone  was  known  to  be  dependent  upon  the  unit  discharge 
rate  and  density  differences  due  to  stratification  and  was  therefore 
considered  to  be  a function  of  the  densimetric  Froude  number.  A form  of 
the  densimetric  Froude  number  used  by  the  Water  Resources  Engineers 
(developed  from  Debler's  criteria)  was  modified  and  used  to  predict  a 
zone  thickness  that  correlated  well  with  thicknesses  observed  in  the 
model.  The  interflow  thickness,  D (ft),  can  be  evaluated  as 


where 

Q = volumetric  flow  rate  of  the  pumpback  current  ( includes  en- 
trainment), cfs 

W = average  reservoir  width  at  the  elevation  'f  the  current  in 
the  vicinity  of  the  dam,  ft 

Ap  = density  difference  of  the  epilimnion  and  hypolimnion,  g/cc 
= average  density  of  the  pumpback  current,  g/cc 

Even  with  unsteady  operation,  the  thickness  of  the  pumpback  current  re- 
mained constant  in  the  vicinity  of  the  dam  for  a particular  stratifica- 
tion, pvunpback  rate,  and  proposed  operational  mode. 


1+9.  Density  profiles  taken  at  the  end  of  each  week  of  prototype 
simulation  revealed  that  the  amount  of  mixing  which  occurred  within  the 
zone  of  interflow  was  proportional  to  time.  This  phenomenon  is  demon- 
strated by  the  temperature  profiles  shown  in  Figure  12  which  are  typical 
of  the  data  obtained  from  the  model.  Model  data  of  this  type  were  used 
to  develop  the  mixing  procedure  and  mixing  coefficients  that  were  used 
in  the  mathematical  model.  A complete  description  of  the  mixing  tech- 
nique is  provided  in  the  mathematical  model  description  (Part  IV). 

temperature 
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PART  IV:  MATHEMATICAL  MODEL  DESCRIPTION 


50.  The  internal  heat  budget  in  Dickey  Lake  and  the  release 
temperatures  downstream  of  Lincoln  School  Lake  were  predicted  using  a 
numerical  simulation  model.  The  model  used  in  conjunction  with  this 
investigation  was  developed  at  WES  to  uniquely  describe  the  thermal  char- 
acteristics a.nd  hydrodynamic  patterns  of  the  Dickey-Lincoln  School  Lakes 

project.  During  the  development  of  the  model,  considerable  insight  was 

0 

gained  from  previous  research  efforts  of  Clay  and  Fruh,  Edinger  and 
Geyer,^  Dake  and  Harleman,^^  Bohan  and  Grace, ^ and  others. 

Fundamental  Assumptions 

51.  The  mathematical  model  provides  a procedure  for  examining 
the  balance  of  thermal  energy  imposed  on  an  impoundment.  This  energy 
balance  and  lake  hydrodynamic  phenomena  are  used  to  map  vertical  pro- 
files of  temperature  and  release  temperatures  in  the  time  domain.  The 
model  includes  computational  methods  for  simulating  heat  transfer  at 

the  air-water  interface;  advective  heat  due  to  inflow,  outflow,  and  pump- 
back  processes;  and  the  internal  dispersion  of  thermal  energy.  ' The 
model  is  one-dimensional  based  on  the  division  of  the  impoundment  into 
discrete  horizontal  layers  of  uniform  thickness.  Assumptions  include 
the  following: 

a.  Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  The  water  in  each  discrete  layer  is  isotropic  and  physi- 
cally homogeneous. 

c_.  Internal  advection  and  heat  transfer  occur  only  in  the 
vertical  direction. 

d.  External  advection  (inflow,  outflow,  and  pmpback)  occurs 
as  a uniform  distribution  within  each  layer. 

£.  Interna],  dispersion  (between  layers)  of  thermal  energy  is 
accomplished  by  a diffusion  mechanism  which  combines  the 
effects  of  molecular  diffusion,  turbulent  diffusion,  and 
thermal  convection. 

52.  The  surface  heat  exchange,  internal  mixing,  and  advection 


36 


processes  are  simulated  separately  and  their  effects  are  introduced 
sequentially  at  specified  time  intervals.  A simplified  flow  chart  of 
the  numerical  simulation  procedure  is  presented  in  Figure  13. 


Siirface  Heat  Exchange 

53.  The  mathematical  model  employs  an  approach  to  the  evaluation 

9 

of  surface  heat  transfer  developed  by  Edinger  and  Geyer.  This  method 
formulates  equilibrium  temperatures  and  coefficients  of  surface  heat 
exchange.  Equilibrium  temperature  is  defined  as  that  temperature  at 
which  the  net  rate  of  heat  exchange  between  the  water  surface  and  the 
atmosphere  is  zero.  The  coefficient  of  surface  heat  exchange  is  the 
rate  at  which  the  heat  transfer  process  will  occur.  The  equation  de- 
scribing this  relationship  is: 

H = K(E  - 0)  (13) 

where 

2 

H = net  rate  of  surface  heat  transfer,  Btu/ft  /day 

2 

K = coefficient  of  surface  heat  exchange,  Btu/ft  /day/°F 
E = equilibrium  temperature,  °F 
6 = siirface  temperature,  °F 

The  compiitation  of  equilibrium  temperature  and  heat  exchange  coefficient 
is  based  solely  on  meteorological  data  as  outlined  by  Edinger,  Dutt- 
weiler,  and  Geyer. 

5^.  The  net  heat  exchange  at  the  surface  is  composed  of  seven 
heat  exchange  processes: 

Shortwave  solar  radiation 
b.  Reflected  shortwave  radiation 
i c_.  Long-wave  atmospheric  radiation 

Reflected  long-wave  radiation 
e^.  Heat  transfer  due  to  conduction 
£.  Back  radiation  from  the  water  surface 
Heat  loss  due  to  evaporation 

For  every  day  of  meteorological  data,  the  seven  heat  exchange  terms  can 
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Figure  13.  Mathematical  model  flow  chart 
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be  evaluated,  the  net  heat  exchange  computed  and  then  expressed  in  terms 
of  an  equilibrium  temperature  and  an  exchange  coefficient. 

55*  All  of  the  surface  heat  exchange  processes,  with  the  excep- 
tion of  shortwave  radiation,  affect  only  approximately  the  top  few  feet 
of  the  lake.  Shortwave  radiation  penetrates  the  water  surface  and  in- 
creases the  temperature  at  greater  depths.  Based  on  laboratory  investi- 
gations, Bake  and  Harleman^*^  have  suggested  an  exponential  decay  with 
depth  for  describing  the  heat  flux  due  to  shortwave  penetration. 

56.  The  siirface  heat  exchange  concepts  are  implemented  in  the 
model  by  the  exponential  penetration  of  a percentage  of  the  incoming 
shortwave  radiation  and  the  placement  of  the  effect  of  all  other  sources 
of  surface  heat  exchange  into  the  surface  layer.  This  proced\ire  can  be 
expressed  mathematically  by  the  following  two  equations: 

H = K(E  - e)  - (1  - 6)S  (li*) 

s 

-Xz . 

H.  = (1  - 6)Se  ^ (15) 

where 

2 

= heat  transfer  rate  into  or  out  of  surface  layer,  Btu/ft  /day 

B = shortwave  radiation  absorbed  in  the  surface  layer,  percent 

2 

S = total  incoming  shortwave  radiation  rate,  Btu/ft  /day 

2 

= rate  of  heat  absorbed  in  layer  i , Btu/ft  /day 
e = natural  logarithmic  base  (2.7183) 

X = heat  absorption  coefficient,  ft  ^ 

= depth  below  surface,  ft 

57-  Equations  lU  and  15  are  applied  to  the  Dickey  Lake  temperature 
profile  once  during  each  one-day  time  step.  A shorter  time  step  was 
evaluated  as  discussed  in  Appendix  A but  was  not  selected  for  use.  The 
net  heat  exchange  rate  into  each  layer  is  computed  and  converted  to  a 
temperature  change.  The  temperature  changes  are  used  to  determine  an  up- 
dated temperature  profile  for  Dickey  Lake.  Numerically,  Lincoln  School 
Lake  is  assumed  to  remain  thoroughly  mixed.  Thus,  surface  heat  exchange 
in  Lincoln  School  Lake  is  effected  by  converting  the  total  heat  exchange 
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rate  as  indicated  in  Equation  13  to  a temperature  change.  The  tempera- 
ture change  is  then  added  to  the  existing  Lincoln  School  Lake  temperature 
once  during  each  one-day  time  step. 

Inflow 

58.  The  process  of  inflow  into  a lake  is  simulated  by  the  place- 
ment of  inflow  quantity  and  quality  at  that  layer  in  which  the  density 
of  the  lake  corresponds  most  nearly  to  the  density  of  the  inflow.  Re- 
search efforts  and  physical  model  studies  at  WES  have  indicated  the 
existence  of  entrainment-induced  density  currents  that  flow  upstream 
along  the  surface  into  the  turbulent  mixing  zone  caused  by  the  inflow. 
Entrainment  is  implemented  in  the  model  by  augmenting  the  inflow 
quantity  with  a volume  from  the  surface  layer.  Characteristics  of  in- 
flow and  entrained  flow  are  averaged;  and  mixed  values  of  density, 
temperature,  and  other  water-quality  parameters  are  determined.  The 
mixed  density  is  used  to  determine  placement  of  the  total  quantity  and 
mixed  quality.  Simulation  of  the  inflow  process  displaces  upward  a 
volume  equal  to  the  total  inflow  quantity.  This  upward  displacement  is 
reflected  in  the  model  by  an  increase  in  the  water-surface  elevation. 

A corresponding  decrease  in  water-surface  elevation  occurs  as  a result 
of  the  outflow  simulation  process. 

59-  The  volume  of  the  entrained  current  is  generally  expressed 
as  a percentage  of  the  inflow  quantity.  Prior  flume  studies  at  WES 
indicated  that  this  percentage  ranges  from  5 to  nearly  250.  The  per- 
centage is  thought  to  be  a function  of  the  entering  velocity,  depth, 
density,  and  density  within  the  lake,  and  is  therefore  a function  of  the 
densimetric  Froude  nunber.  However,  analytical  relationships  have  not 
been  completed. 

60.  As  indicated  previously  in  the  discussion  of  the  Dickey- 
Lincoln  School  Lakes  distorted  model,  the  Inflow  current  in  Dickey  Lake 
has  no  appreciable  effect  on  the  temperature  profile  in  the  vicinity  of 
Dickey  Dam.  The  same  conclusion  was  reached  numerically  by  simulation 
with  various  inflow  temperatures.  Simulations  with  inflow  temperatures 
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(a)  as  provided  by  NED,  (b)  with  a constant  value  of  0°C  throughout  the 
simulation  period,  or  (c)  with  a constant  value  of  100°C  produced  iden- 
tical temperature  profiles  within  the  accuracy  of  the  model  output.  Thus, 
it  was  shown  that  due  to  the  large  volume  of  Dickey  Lake,  the  inflow  pro- 
cess has  a minimal  effect  on  the  heat  budget  of  the  lake  at  the  dam. 
However,  the  inflow  process  must  be  retained  in  the  model  in  order  to 
maintain  the  water  budget.  The  minimal  effect  of  inflow  on  the  heat 
budget  allowed  the  flow  quantity  from  many  inflow  points  to  be  lumped 
together  and  treated  as  a single  inflow  point  with  total  flow  quantity 
and  flow-weighted  average  temperature.  Lumping  of  inflows  decreases 
computer  storage  requirements  and  run  time. 


Internal  Dispersion 


6l.  The  internal  dispersion  process  is  represented  by  an  internal 
mixing  scheme  based  on  a simple  diffusion  analogy.  Internal  mixing 
transfers  heat  and  other  water-quality  constituents  between  adjacent 
layers.  The  magnitude  of  the  transfer  between  two  layers  is  a percent- 
age of  the  total  transfer  required  to  completely  mix  the  two  layers. 

This  percentage  is  a mixing  coefficient  which  is  defined  for  every  layer. 
Data  input  includes  values  of  the  mixing  coefficient  at  the  top  and  at 
the  bottom  of  the  lake.  An  exponential  fit  between  the  two  extreme 
values  is  used  to  determine  the  appropriate  coefficient  at  each  layer. 


Outflow 

62.  The  outflow  component  of  the  model  incorporates  the  selective 
withdrawal  predictive  techniques  developed  at  WES^  for  submerged  weir 
flow.  However,  the  selective  withdrawal  description  was  modified  to  in- 
clude characteristics  unique  to  the  Dickey  Lake  intake  structure  (para- 
graphs 17-21.  Transcendental  equations  defining  the  location  of  the  zero- 
velocity  limits  are  solved  with  a half-interval  search  method.  The  loca- 
tion of  the  zero-velocity  limits  is  functionally  dependent  on  config\ira- 
tion  of  the  withdrawal  device,  release  flow  rate,  and  density  structure 
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within  the  lake.  With  knowledge  of  the  limits  of  the  withdrawal  zone, 
the  velocity  profile  due  to  outflow  can  be  determined.  The  flow  from 
each  layer  is  then  computed  as  the  product  of  the  velocity  in  the  layer 
and  the  width  and  thickness  of  the  layer.  A flow-weighted  average  is 
applied  to  the  temperature  profile  to  determine  the  value  of  the  release 
temperature  for  each  time  step.  For  the  constant  discharge  rates  pro- 
posed during  periods  of  generation,  the  physical  models  indicated  that 
the  withdrawal  zone  was  fully  established  almost  immediately  after  initi- 
ation of  withdrawal.  Additionally,  release  densities  were  constant 
throughout  a generation  period.  Therefore,  in  the  mathematical  model  the 
release  temperature  is  assumed  constant  throughout  a generation  period. 


Operation  Schedules 


63.  Because  the  generation  and  pumpback  flows  change  markedly 
within  the  period  of  a day  (Table  l)  and  the  rates  of  these  flows  affect 
withdrawal  and  pumpback  characteristics,  it  was  not  adequate  to  use 
daily  average  flows  for  generation  and  pumpback  operations.  Daily  aver- 
age flow  routings  are  quite  adequate  for  less  dynamic  reservoirs , such 
as  water  supply  and  flood  control  projects.  For  this  model,  the 
scheduled  generation  and  pumpback  flow  rates  and  durations  are  used 
rather  than  daily  averages.  The  Julian  day  number,  the  modes  of  opera- 
tion (generation  or  pumpback),  flow  rates,  durations  of  flow,  and  eleva- 
tions of  the  withdrawal  or  pumpback  device  are  input  for  each  day  that 
a change  in  the  operation  schedule  occurs.  The  flow  rates  are  used  in 
computing  withdrawal  and  pumpback  characteristics,  while  the  duration 
and  rate  of  flow  are  used  to  budget  generation  and  pumpback  volumes. 


Lincoln  School  Lake 

6U.  Lincoln  School  Lake  is  modeled  numerically  by  maintaining 
heat  and  water  budgets  in  the  lake  and  is  assimied  to  remain  fully 
mixed  throughout  the  year.  This  assumption  is  supported  by  physical 
model  observation,  and  it  allows  further  simplifying  assumptions  to 
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be  made  with  respect  to  advection  and  heat  exchange  in  the  lake.  Once 
during  each  day  of  simulation  the  volume  of  water  in  the  lake  is  ad- 
justed to  sequentially  account  for  (a)  Allagash  River  volume  which  is 
not  pumped  back,  (b)  Dickey  Lake  generation  volume  which  is  not  pumped 
back,  (c)  Lincoln  School  Lake  volume  which  is  pumped  back,  and  (d)  Lin- 
coln School  Lake  volume  which  is  released  downstream.  The  net  contri- 
butions of  volume  from  the  Allagash  River  and  Dickey  Lake  at  the  respec- 
tive temperatures  are  added  to  the  Lincoln  School  Lake  volume,  producing 
a mixed  value  of  temperature.  The  daily  heat  budgeting  is  completed  by 
applying  Equation  13  to  the  final  surface  area  of  Lincoln  School  Lake  to 
account  for  surface  heat  exchange  at  the  air-water  interface.  The 
schematic  shown  in  Figure  ih  describes  the  budgeting  of  heat  and  mass  in 
Lincoln  School  Lake. 
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Figure  ih.  Simulation  of  reregulation 
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Pumpback  Quality 

65.  The  quality  (temperature)  of  the  pumped  flow  can  be  in- 
fluenced by  three  water  sources;  the  Allagash  River,  releases  from 
Dickey  Lake  just  prior  to  pumpback,  and  Lincoln  School  Lake.  As  pre- 
viously discussed  (paragraph  UU),  the  Allagash  River  was  found  to  con- 
tribute to  the  quality  of  pumped  flow.  The  amount  of  contribution  is 
described  by  Equation  11.  The  remainder  of  the  water  pumped  back  can 
consist  of  either  Dickey  Lake  release  or  Lincoln  School  Lake  water. 

If  generation  occurs  before  pumpback  on  the  same  day,  then  the  remainder 
will  consist  of  Dickey  Lake  release  water;  otherwise,  Lincoln  School 
Lake  water  constitutes  the  remainder.  The  temperature  of  the  pumped 
flow  is  determined  from  a volume-weighted  average  of  the  constituents. 

Pumpback  Mixing  Technique 

66.  As  previously  indicated  in  the  description  of  the  physical 
model  results,  the  pumpback  current  is  a two-dimensional  phenomenon 
(Figure  ll)  that  causes  longitudinal  variations  in  the  density  structure 
of  Dickey  Lake.  Different  thermal  profiles  existed  in  the  vicinity  of 
the  dam  than  those  observed  farther  upstream  in  the  lake.  As  the  math- 
ematical model  is  one-dimensional,  the  decision  was  made  to  simulate  the 
thermal  profile  immediately  upstream  of  the  dam,  since  this  profile 
controls  the  withdrawal  and  pumpback  characteristics  within  Dickey  Lake. 

67.  In  the  mathematical  model,  the  pumpback  current  is  placed  in 
the  layer  whose  density  most  nearly  corresponds  to  the  density  of  the 
pumpback  current.  In  the  vicinity  of  the  dam,  the  pumpback  current 
actually  affects  a zone  above  and  below  this  layer  (paragraph  U8).  To 
represent  this  effect  and  still  maintain  the  heat  budget  of  the  lake, 
the  layers  within  the  zone  are  mixed.  The  mixing  technique  employed  is 
based  on  the  concept  of  a portion  of  each  of  the  layers  in  the  pumpback 
zone  being  removed,  mixed  together,  returning  the  mixed  water  to  the 
layers  in  the  zone,  and  then  mixing  within  each  layer.  The  portion  of 
water  removed  and  mixed  from  each  layer  is  computed  by  multiplying  a 
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mixing  coefficient  (a  number  between  0 and  l)  times  the  volume  of  the 
layer.  The  mixing  coefficients  were  determined  using  data  from  the  dis- 
torted physical  model  (paragraph  h9) • Density  profiles  were  measixred 
before  and  after  pumpback.  By  analyzing  the  difference  in  these  pro- 
files, the  mixing  coefficients  necessary  to  apply  to  the  initial  profile 
to  yield  the  final  profile  were  computed. 

68.  Mixing  coefficients  were  computed  for  each  of  the  pumpback 
operation  tests  conducted  in  the  distorted  model.  Analysis  of  the 
mixing  coefficients  revealed  that  the  coefficients  could  be  considered 
symmetrical  about  the  pvmpback  current  inflow  layer.  Additionally,  a 
constant  magnitude  of  the  coefficients  could  be  used  throughout  the 
stratification  cycle.  For  a given  pumping  rate,  the  mixing  coefficients 
could  be  assumed  equal  for  pumpback  over  or  under  the  weir.  For  three 
reversible  units  pumping  20,000  cfs,  the  coefficients  were  found  to  be 
approximately  twice  the  magnitude  of  those  determined  for  pumpback  with 
a single  reversible  turbine  (7,000  cfs).  This  difference  was  attributed 
to  the  greater  velocity  of  the  current  with  pumped  flow  of  20,000  cfs. 

An  equation  was  developed  to  predict  the  mixing  coefficient  in  a layer 
as  a function  of  the  distance  of  that  layer  from  the  pumpback  current 
inflow  layer: 

-BX. 

n.  = Ae  ^ (l6) 

where 

= mixing  coefficient 
i = layer 

A = 0.067  for  3 reversible  units 

= 0.033  for  1 reversible  unit 

B = 0.118 

= distance  from  the  ptunpback  current  inflow  layer  to  layer,  i 

69.  To  represent  the  pumpback  mixing  technique  within  the  mathe- 
matical model,  the  following  sequence  of  computational  steps  is  em- 
ployed. The  depth  of  the  zone  affected  by  the  pumpback  current  is 
computed  by  Equation  12,  and  this  zone  is  distributed  about  the'layer 
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of  piimpback  placement.  The  mixing  coefficient  for  each  layer  in  the 
zone  is  computed  by  Equation  l6,  the  amount  of  each  layer  to  be  mixed 
is  computed  as  the  product  of  the  mixing  coefficient  and  the  layer 
volume,  and  the  portions  from  the  layers  are  mixed  together.  The  mixed 
water  is  returned  to  the  layers  and  each  layer  is  individually  mixed. 
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PART  V:  DISCUSSION 

TO.  A one-dimensional,  numerical  simulation  model  was  developed 
for  prediction  of  temperature  within  and  downstream  of  Dickey-Lincoln 
School  Lakes.  Physical  hydraulic  models  were  used  to  obtain  an  improved 
understanding  and  description  of  the  hydrodynamic  response  of  the  lakes 
to  pumped-storage  hydropower  operations.  Descriptions  from  the  physical 
models  were  incorporated  into  the  mathematical  model  to  enhance  the  pre- 
dictive capabilities.  Sensitivity  analyses  were  conducted  to  determine 
the  effects  on  temperature  of  the  various  algorithms  developed  from  the 
physical  model  studies.  Discussions  of  the  sensitivity  analyses  are  pre- 
sented in  the  following  paragraphs.  It  should  be  noted  that  the  results 
of  these  sensitivity  analyses  are  unique  to  the  Dickey-Lincoln  School 
Lakes  project. 

Tl.  Information  obtained  from  the  Dickey  Lake  intake  model 
(Part  II)  was  used  to  modify  the  mathematical  description  of  the  selec- 
tive withdrawal  characteristics  (paragraphs  17-21).  The  effect  that 
these  modifications  have  on  the  predicted  release  temperatures  from 
Dickey  and  Lincoln  School  Lakes  is  presented  in  Plates  5a  and  5b,  re- 
spectively. The  dashed-line  circle  represents  a base  condition  of  re- 
lease temperature  computed  without  the  use  of  the  selective  withdrawal 
modifications.  The  solid  line  plots  the  deviation  in  temperature  from 
the  base  condition  as  a result  of  the  selective  withdrawal  modifications. 
The  base  and  perturbed  conditions  are  plotted  over  the  cycle  of  a year 
(360°  = 1 year)  starting  on  1 January  positioned  at  the  scaled  axis  and 
rotating  in  a counterclockwise  direction.  Hydrologic  and  meteorologic 
data  for  1953,  furnished  by  NED,  and  operational  conditions  for  six 
turbines  (generation  flow  60,000  cfs  and  pumpback  flow  20,000  cfs.  Ta- 
ble 1)  were  used  to  conduct  all  the  sensitivity  analyses.  As  shown  by 
Plate  5a,  the  selective  withdrawal  modifications  of  the  mathematical 
model  had  the  effect  of  decreasing  the  Dickey  Lake  release  temperature 
by  a maximum  of  about  3°C.  The  reregulation  pool  attenuated  this  effect 
to  a maximum  of  approximately  2°C  (Plate  5b).  Although  these  modifica- 
tions did  not  create  drastic  changes  in  the  release  temperature,  the 
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changes  are  significant  and  cannot  be  overlooked. 


72.  The  ability  of  the  reregulating  pool  to  warm  and  cool  the  re- 
leases from  Dickey  Lake  is  demonstrated  by  Plate  6.  The  surface  heat 
exchange  expected  in  Lincoln  School  Lake  tends  to  drive  the  water  tem- 
perature toward  equilibrium  with  the  atmosphere.  During  the  strong 
stratification  period,  the  releases  from  Lincoln  School  Lake  should  be 
as  much  as  5°C  warmer  than  the  Dickey  Lake  release.  Other  factors 
contributing  to  the  heat  exchange  characteristics  of  the  reregulation 
pool  are  the  influence  of  the  Allagash  River,  the  contribution  of  the 
Dickey  release,  the  amount  of  release  and  pumpback  from  Lincoln  School, 
and  the  fully  mixed  conditions  of  the  lake  (Figure  lb).  Descriptions  of 
these  water  budgeting  processes  in  Lincoln  School  Lake  were  developed 
from  hydrodynamic  model  results  (Part  III)  for  use  in  the  mathematical 
model.  For  example,  the  hydrodynamic  model  was  used  to  determine  the 
contribution  of  the  Allagash  River  to  the  pumped  flow  from  Lincoln 
School  Lake  (paragraphs  39  and  bb).  The  Allagash  River  water  that  is 
not  pumped  back  contributes  to  the  heat  budget  of  Lincoln  School  Lake. 
The  difference  in  the  Lincoln  School  release  temperature  while  allowing 
the  pumpback  of  Allagash  River  water  (using  Equation  ll)  compared  with 
not  allowing  any  pumpback  of  Allagash  water  is  plotted  in  Plate  T.  This 
plot  demonstrates  the  importance  of  providing  a reasonable  mathematical 
description  of  the  advected  heat  into  or  out  of  Lincoln  School  Lake. 

73-  Physical  understanding  of  the  pumped  flow  characteristics  was 
obtained  from  both  the  undistorted-  and  distorted-scale  model  results. 
Two  basic  features  of  the  pumped  flow  are  entrainment  of  Dickey  Lake 
water  into  the  pumpback  jet  and  mixing  within  Dickey  Lake  created  by  the 
pumpback  current.  A description  of  these  features  was  considered  nec- 
essary to  provide  reasonable  predictions  of  thermal  profiles  within 
Dickey  Lake.  The  entrainment  and  pumpback  mixing  algorithms  of  the 
mathematical  model  did  affect  the  development  of  thermal  profiles  in 
Dickey  Lake.  However,  the  effect  of  the  entrainment  phenomena  (de- 
scribed in  paragraphs  25  and  26)  on  the  release  temperature  from  Dickey 
Lake  was  very  small  (Plate  8a)  and  the  effect  on  the  release  from 
Lincoln  School  Lake  was  practically  undetectable  (Plate  8b).  With  a 
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lake  as  large  as  Dickey,  entrainment  had  a small  effect  on  the  thermal 
profile  which  resulted  in  an  almost  unnoticeable  difference  in  the  re- 
lease temperature.  The  pumphack  mixing  algorithm  discussed  in  Part  IV 
(paragraphs  66-69)  had  an  apparent  warming  effect  (Plate  9a)  on  the 
Dickey  release  compared  with  the  base  condition  of  no  mixing  of  the 
pumpback  current  within  the  Dickey  thermal  profile.  This  warming  trend 
was  damped  in  Lincoln  School  Lake  as  shown  by  Plate  9b.  The  pumpback 
current  mixing  technique  resulted  in  a deepening  of  the  metalimnion, 
thus  providing  a larger  tone  for  the  release  of  warm  water.  The  spikes 
on  the  plots  of  Plates  8 and  9 are  due  to  a mathematical  anomaly  re- 
sulting from  a slight  lead  in  the  development  of  the  thermal  profile 
for  the  perturbed  condition  as  compared  with  the  base  condition. 

7^.  Simulations  were  conducted  with  pumpback  over  the  gates  and 
under  the  gates  to  evaluate  the  effect  of  pumpback  location  on  release 
temperatures.  Although  some  variation  in  the  temperature  profiles  for 
the  two  conditions  was  observed,  the  differences  in  the  release  tempera- 
tures were  negligible  as  shown  by  Plate  10.  Therefore,  these  results 
indicate  that  there  is  no  advantage  to  pumping  back  lower  in  the  pool 
versus  pumping  back  high  (over  the  weir). 

75.  Definition  of  the  selective  withdrawal  characteristics  of 
Dickey  Lake  had  by  far  the  greatest  impact  on  release  temperatures. 
Variations  in  the  pumpback  characteristics  only  slightly  affected  re- 
lease temperatirres . This  was  attributed  to  the  large  volume  of  Dickey 
Lake  and  large  zone  of  withdrawal  that  existed  with  the  generation 
flows.  The  release  temperature  is  a flow-weighted  function  of  the  tem- 
peratures within  the  zone  of  withdrawal;  however,  release  temperature 
is  less  sensitive  to  variations  in  the  temperature  profile  for  a large 
withdrawal  zone.  For  different  release  conditions  causing  a smaller 
withdrawal  zone,  variations  in  the  temperature  profile  due  to  pumpback 
could  cause  more  significant  differences  in  release  temperatures. 

76.  Simulations  were  conducted  with  the  effects  of  the  restric- 
tive topography  on  selective  withdrawal  removed  (paragraphs  16,  I8, 
and  19)  to  determine  the  variation  in  release  temperature  from  the  base 
condition  which  includes  the  topographical  effect.  The  perturbation 
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plot  of  Plate  11a  demonstrates  that  the  release  temperatures  from  Dickey 
Lake  are  increased  1/2  to  2°C  during  the  peak  stratification  season. 

This  increase  is  damped  in  Lincoln  School  Lake  to  an  increase  of  1/2  to 
1°C  (Plate  11b).  Comparison  of  Plates  5a  and  11a  shows  that  the  topo- 
graphical effect  accounts  for  about  half  of  the  total  variation  from 
the  condition  of  no  selective  withdrawal  modifications.  With  excessive 

topographical  contraction,  selective  withdrawal  can  be  completely  r re- 
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vented  as  shown  by  Kao.  As  mentioned  previously,  the  topograf. r..y 
partially  influenced  selective  withdrawal  from  Dickey  Lake,  ani 
tive  withdrawal  was  not  eliminated. 

TT.  Results  of  this  model  study  and  sensitivity  analyses  are 
unique  to  the  Dickey-Lincoln  School  Lakes  project  and  should  not  be 
interpreted  to  apply  to  other  pumped-storage  projects. 
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Operational  Modes  Tested 


Lincoln  School  base  generation  flow  was  1000  cfs  for  times  not  noted  above.  The  duration  eind 
frequency  of  releases  from  Lincoln  School  Lake  were  dependent  upon  the  magnitude  of  the  Allagash 
inflow  that  was  being  simulated. 
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APPENDIX  A:  TIME  STEP  FOR  HEAT  EXCHANGE  COMPUTATIONS 


1.  The  purpose  of  this  analysis  was  to  determine  the  appropriate 
time  step  for  applying  the  surface  heat  exchange  in  the  numerical  simu- 
lation model.  Conventionally,  one-dimensional  reservoir  thermal  simu- 

13* 

lation  models  employ  a daily  time  step.  Harleman  and  Hurley  compared 
three  hourly  and  dally  surface  heat  exchange  inputs  in  a one-dimensional 
model.  They  found  variations  in  surface  temperatures  of  2-3°C  over  a 
day.  Additionally,  they  noted  that  diurnal  temperature  fluctuations 
affected  primarily  the  upper  2-k  m of  the  lake.  Their  results  showed 
that  the  sirrface  temperatures  and  thermal  profiles  based  on  thr’ee  hourly 
inputs  fluctuated  about  the  values  predicted  with  average  daily  input. 
Since  the  mathematical  model  for  Dickey-Lincoln  School  Lakes  incorpo- 
rates the  simulation  of  generation  and  pumpback  with  a time  step  less 
than  one  day,  the  need  for  a similar  time  step  for  surface  heat  exchange 
input  was  evaluated. 

2.  Numerical  simulations  of  Dickey  Lake  were  made  using  an  hourly 
and  a daily  time  step  for  surface  heat  exchange  input.  Plate  AI  shows 
daily  values  and  mean  hourly  values  i2k-hr  period)  of  Dickey  Lake  re- 
lease temperatures.  Plate  A2  compares  typical  thermal  profiles  for  the 
daily  time  step  with  selected  hourly  profiles.  Plate  A3  shows,  for  one 
simulation  day,  the  fluctuations  of  the  hourly  surface  temperatures 
about  the  daily  value.  Plate  kk  shows,  for  a stratification  period, 
the  fluctuations  of  the  mean  hourly  release  temperatures  compared  with 

the  daily  values.  Analysis  of  these  results  yielded  conclusions  simi- 

13 

lar  to  those  of  Harleman  and  Hurley.  The  effects  of  the  hourly  time 
step  were  small  and  restricted  to  the  upper  portion  of  the  profile. 

The  effect  on  the  release  temperature  from  Dickey  Lake  was  minimal. 

3.  Based  on  the  results  of  the  simulation,  a daily  time  step  was 
chosen  for  surface  heat  exchange  input  for  the  Dickey-Lincoln  School 
numerical  model.  Since  this  model  computes  the  effect  of  individual 
processes  on  the  total  heat  budget  of  ^he  lake  and  then  applies  these 


* See  references  at  end  of  main  text. 
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effects  sequentially,  the  model  can  easily  handle  different  time  steps 
for  each  process.  This  means  that  a daily  time  step  for  surface  heat 
exchange  input  and  a shorter  time  step  for  generation  or  pumpback  can 
be  used.  Additionally,  the  use  of  a daily  time  step  versus  a shorter 
time  step  for  surface  heat  exchange  input  results  in  a computational 
time  savings  and  a reduction  in  input  data  requirements.  If  surface 
heat  exchange  input  were  to  be  required  for  a time  step  less  than  daily, 
that  input  (equilibrium  temperatures  and  exchange  coefficients)  would 
have  to  be  derived  from  daily  values.  This  derivation  would  be  accom- 


plished by  curve-fitting  an  assumed  functional  relationship  to  the  daily 
ll 

values.  The  computed  hourly  input  may  or  may  not  be  representative  of 


the  actual  hourly  fluctuations.  Finally,  these  functional  relationships 
ensure  that  the  hourly  values  fluctuate  symmetrically  about  the  daily 
average  value. 
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APPENDIX  B;  NOTATION 
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A 0.06T  for  3 reversible  units;  0.033  for  1 reversible  unit 
b Weir  length,  ft 

B 0,ll8 

D Interflow  thickness  at  a particular  station,  ft 
e Natural  logarithmic  base  (2.7183) 

E Equilibrium  temperature,  °F 

2 

g Acceleration  due  to  gravity,  32.2  ft/sec 

2 

H Net  rate  of  surface  heat  transfer,  Btu/ft  /day 

2 

H.  Rate  of  heat  absorbed  in  layer  i , Btu/ft  /day 

^ 2 

Heat  transfer  rate  into  or  out  of  surface  layer,  Btu/ft  /day 

H Head  on  the  weir  for  free  flow  or  depth  of  flow  over  the 
weir  for  submerged  flow,  ft 

i Layer 

2 

K Coefficient  of  surface  heat  exchange,  Btu/ft  /day/°F 

p Exponent  in  the  selective  withdrawal  dimensionless  velocity 
distribution  equation 

Q Discharge  over  weir,  cfs 

Volume  flow  rate  of  the  pumpback  current,  cfs 

Q Volume  rate  of  the  entrained  flow,  cfs 

Volume  flow  rate  of  water  pumped,  cfs 

Q Allagash  River  inflow,  1000  acre-ft/day 

A 2 

S Total  incoming  shortwave  radiation  rate,  Btu/ft  /day 

v^  Local  velocity  in  the  zone  of  withdrawal  at  a distance  y^ 

below  the  elevation  of  the  maximum  velocity  V , fps"^ 

max 

Vg  Local  velocity  in  the  zone  of  withdrawal  at  a distance 

above  the  elevation  of  the  maximum  velocity  V , fps 

•'  max  ^ 

V Maximum  velocity  in  the  zone  of  withdrawal,  fps 

iTid.x 

V Average  velocity  over  the  weir,  fps 
w 

W Average  reservoir  width  at  the  elevation  of  the  pumpback 
current  in  the  vicinity  of  the  dam,  ft 

X.  Vertical  distance  from  the  pumpback  current  inflow  layer  to 
layer  i , ft 

y^  Vertical  distance  from  the  elevation  of  the  maximum  velocity 
V to  that  of  the  corresponding  local  velocity  v , ft 
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Vertical  distance  from  the  elevation  of  the  maximum  velocity 
V to  that  of  the  corresponding  local  velocity  v , ft 

iT13JC  c. 

Vertical  distance  from  the  elevation  of  the  maximum  velocity 
to  the  lover  limit  of  withdrawal,  ft 
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At 

Ap 
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Ap 


Im 


Ap 
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Vertical  distance  from  the  elevation  of  the  maximum  velocity 
to  the  upper  limit  of  withdrawal  or  the  free  surface  for 
weir  flow,  ft 

Modified  vertical  distance  from  the  elevation  of  the  maximum 
velocity  to  the  lower  limit  of  withdrawal,  ft 

Depth  below  surface,  ft 

Vertical  distance  from  the  elevation  of  the  weir  crest  to  the 
lower  limit  of  the  zone  of  withdrawal,  ft 

Shortwave  radiation  absorbed  in  the  surface  layer,  percent 

Pumping  period 

Density  difference  of  the  epilimnion  and  hypolimnion,  g/cc 

Density  difference  of  fluid  between  the  elevations  of  the 
maximum  velocity  and  the  corresponding  local  velocity  v^  , 
g/cc 

Density  difference  of  fluid  between  the  elevations  of  the 
maximum  velocity  and  the  lower  limit  of  the  zone  of  with- 
drawal, g/cc 

Density  difference  of  fluid  between  the  elevations  of  the  weir 
crest  and  the  lower  limit  of  the  zone  of  withdrawal,  g/cc 

Entrainment  coefficient 

Mixing  coefficient 

Surface  temperature,  °F 

Heat  absorption  coefficient,  ft  ^ 

Average  density  of  the  pumpback  current,  g/cc 

Average  density  of  the  entrained  flow,  g/cc 

Density  of  the  pumped  water,  g/cc 

Density  of  fluid  at  the  elevation  of  the  weir  crest,  g/cc 
Percent  of  Allagash  River  water  to  the  total  pumped  flow 
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